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The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi-step process which involves data from 
three different telescopes: the 1.0-m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which 
was used to establish the SDSS standard star network); the SDSS 0.5-m Photometric Telescope (PT) at the Apache Point 
Observatory (APO), New Mexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and 
the SDSS 2.5-m telescope at APO (which obtains the imaging data for the SDSS proper). 

In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data 
from the single-CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0-m and the PT). We 
also describe transformation equations that convert photometry on the USNO- 1.0m u'g'r'i'z' system to photometry the 
SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the 
semi-automated PT factory, which runs MTPIPE in the day-to-day standard SDSS operations at Fermilab. Finally, we 
discuss the use of MTPIPE in current SDSS-related projects, including the Southern u'g'r'i'z' Standard Star project, the 
u'g'r'i'z' Open Star Clusters project, and the SDSS extension (SDSS-II). 
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1 Introduction 

The Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Stoughton et al. 2002a; Abazajian et al. 2003, 2004, 2005; 
Adelman-McCarthy et al. 2006) is a modern, CCD-based 
optical imaging and spectroscopic survey of the Northern 
Galactic Cap. Both imaging and spectroscopy are performed^ 
using a 2.5m f/5 Ritchey-Chretien telescope (Gunn et al. 
2006). The imaging camera (Gunn et al. 1998) contains an 
imaging array of thirty 2048x2048 SITe CCDs and scans 
the sky in drift scan mode along great circles in five dif- 
ferent filters (ugriz). A complete scan from the western 
to the eastern borders of the survey area is called a strip. 
Since there are gaps within the CCD mosaic, it requires 
two strips, offset by 93% of a CCD width, to fill in a com- 
plete rectangular area (which is called stripe). The imaging 
portion of the SDSS has a depth of r = 22.2 (95% detec- 
tion repeatablity for point sources; Ivezic et al. 2004). Fol- 
lowup spectrosopy is performed for a variety of targetted 
galaxies (Eisenstein et al. 2001; Strauss et al. 2002), quasars 
(Richards et al. 2002), and stars via a pair of 320-fiber multi- 
object spectrographs. 

The SDSS imaging covers several thousand square de- 
grees of sky, and over this region the SDSS photometric cal- 
ibrations achieve an accuracy of «0.02 mag (2%) rms in g, 
r, and i and «0.03 mag (3%) in u and z (Ivezic et al. 2004; 
Adelman-McCarthy et al. 2006). The photometric calibra- 
tion of the SDSS 2.5m imaging data is a multi-step process 
which involves the data from three different telescopes and 
three different data processing pipelines. The telescopes in 
question are: 

1. the US Naval Observatory (USNO) 1.0-m telescope at 
Flagstaff Station, Arizona, which was used to set up a 
network of 158 primary standard stars for the u'g'r'i'z' 
photometric system in the magnitude range 9 ^ r' Si 14 
(Fig-Cl Fukugita et al. 1996; Smith et al. 2002); 

2. the 0.5-m Photometric Telescope (PT) at Apache Point 
Observatory (APO), New Mexico, which observes a set 
of u'g'r'i'z' primaries over a range of airmasses to de- 
termine the photometric solution for the night (zeropoints| 
& extinctions), and calibrates stars down to r' w 18 in 
transfer fields — or secondary patches — that are placed 
throughout the SDSS survey area (Fig. [2J; and 

3. the SDSS 2.5-m telescope itself (Gunn et al. 2006), also 
located at APO, whose imaging camera (Gunn et al. 
1998) scans over the secondary patches (Fig.|2ji during 
normal course of operations. Note that the imaging cam- 
era saturates under normal operating conditions at r sa 
14, necessitating the use of faint stars in the secondary 
patches (rather than observing the primary standards di- 
rectly). The secondary patches, which are 40' x 40' in 
size, are grouped in sets of four, so that each group spans 
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the full 2.5° width of a survey stripe, and each group 
is spaced at roughly 15° (1 hour) intervals along a sur- 
vey stripe. [The positions of patches along a stripe have 
been optimized to avoid bright (V SS 8) stars; hence the 
slight variability in the spacings between patch groups 
and between patches within patch groups.] To calibrate 
the 2.5m imaging camera photometry, atmospheric ex- 
tinction measurements are taken from PT observations 
for the night, and photometric zeropoints are obtained 
by matching stars in the imaging camera data with those 
from the secondary patches. This method of calibration 
permits the 2.5m telescope to spend its time more effi- 
ciently, in that the onus of obtaining nightly extinction 
measurements and of calibrating primary and secondary 
standard stars is shifted to other, smaller telescopes. 

The three software pipelines are: 

1. the Monitor Telescope Pipeline (MTPIPE) 1 , which was 
used to process the USNO- 1.0m data for setting up the 
u'g'r'i'z' primary standard star network (Smith et al. 
2002), and which is currently used to process the PT 
data for determining the atmospheric extinctions at APO 
each night and for calibrating the secondary patch fields 
to the SDSS ugriz system; 

2. the Photometric Pipeline (PHOTO; Lupton et al. 2001; 
Lupton 2006), which processes the 2.5m imaging cam- 
era data, yielding among its many outputs accurate as- 
trometry (Pier et al. 2003) and ugriz instrumental pho- 
tometry (counts); and 

3. the (New) Final Calibrations Pipeline (NFCALIB; § 4.5.3.| 
of Stoughton et al. 2002a), which matches the stars in 
the PT secondary transfer fields with the 2.5m imag- 
ing camera photometry, calculates the resulting photo- 
metric zeropoints, and applies them to the 2.5m instru- 
mental magnitudes. (To deal with 2.5m imaging scans 
that are observed over a range of airmasses, NFCALIB 
also makes use of the first-order atmospheric extinction 
coefficients measured by the PT on the same night as 

a given imaging scan. Since most 2.5m imaging scans 
only cover a small range in airmass, this is typically a 
small effect.) 

These pipelines were written with the expressed purpose 
of meeting the image processing and calibration needs of 
the SDSS, and are in no way meant to supplant the many ex- 
cellent general-purpose astronomical image processing and 
analysis packages like IRAF 2 (Tody 1986, 1993), MIDAS 3 
(Warmels 1985; Grosb0l 1989), Gypsy 4 (Shostek & Allen 
1980; Allen et al. 1992; van der Hulst et al. 1992), SExtrac- 
tor (Berlin & Arnouts 1996), DAOPHOT (Stetson 1987), 

1 The Monitor Telescope Pipeline is named after the original SDSS 
0.6m calibration telescope at APO, called the Monitor Telescope (or MT), 
which was de-commissioned in Summer 1998 and replaced by the PT, 
which was commissioned in Spring 1999. Overlapping this time frame, 
the USNO-l.Om telescope was used to set up the u'g'r'i'z' standard star 
network, a task originally planned for the old MT. 

2 iraf.noao.edu 

3 |http : / / www . eso . org/pro jects/ esomidas /| 

4 |http : / /www, astro . rug ■ nl/~gipsy/| 
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* Northern + Equatorial Standards 




Fig. 1 Locations of the 158 Smith et al. (2002) primary standards (asterisks) and the 64 fields of the southern extension 
to the primary standard star network (unfilled squares; see §|5J. Circled symbols indicate u'g'r'i'z' standard stars or fields 
for which there is currently Landolt U BVR C I C photometry (Landolt 1973, 1983, 1992) or for which Landolt is currently 
obtaining UBVR C I C photometry (Landolt, in prep.). 




Fig. 2 Locations of the SDSS secondary patches in equatorial coordinates (Aitoff projection): (black) Northern SDSS 
patches, (red) Southern SDSS patches, and (green) SEGUE patches. (SEGUE is part of the Sloan extension, or SDSS-II, 
and is discussed in §|5]) 



and DoPHOT (Schechter et al. 1993), to name but a few. 
Thus, the SDSS imaging and photometric calibration pipelin 
have more in common with such special-purpose software 
packages as the ESO Imaging Survey pipeline (Nonino et 
al. 1999), the INT Wide Field Camera pipeline (Irwin & 
Lewis 2001), the VLT Survey Telescope pipeline (Grado et 
al. 2004), and the Liverpool Telescope Gamma Ray Burst 
pipeline (Guidorzi et al. 2006). 

The topic of the current paper is MTPIPE. In the fol- 
lowing sections, we will discuss in turn MTPIPE's compo- 



nent packages (§ EJ, validation tests of the MTPIPE soft- 
ware (§|5J, the MTPIPE-based semi-automated PT data pro- 
cessing factory (§ |4}, and other SDSS-related projects us- 
ing MTPIPE (§0; in §[6] we comment on future plans for 

MTPIPE. 

In closing, we note that, although the SDSS magnitudes 
and colors are shown without primes (ugriz), the magni- 
tudes of the primary standard stars are shown with primes 
(v! g'r'i' z'). This is no accident. The primed system is de- 
fined in the natural system of the USNO-l.Om, its CCD, 
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full set of frames (u'g'r'i'z 1 ) are present for each target se- 
quence, and creates histograms for each of the bias and flat 
field frames. 

As an example, plotted in Figure|4]are histograms show- 
ing the distribution of pixel values for each of the bias frames| 
observed by the PT on MJD 53501. 5 Note that the PT CCD 
uses two-amplifier readout, as evidenced by the strongly bi- 
modal distribution of pixel values for these bias frames. 

Note that MTPIPE does not use the values of image 
FITS header keywords directly, but makes use of an ASCII 
observing log (called the mdReport file) containing the rel- 
evant information (e.g., RA, DEC, time of observation, ex- 
posure time, filter, ...) for each image. This information is 
usually based upon image FITS header keyword values, but 
can also be user-generated. The use of the mdReport file 
simplifies the process of editing this information by humans 
when the need arises. 

2.2 mtFrames 

The next package, mtFrames, is the primary image pro- 
cessing portion of the software, and is capable of processing 
data from CCDs with 1 -, 2-, or 4-amplifier readout electron- 
ics. 

The first part of mtFrames creates master bias frames, 
flat field frames, and fringe frames. The master biases are 
created by median filtering a set raw bias frames. During 
general operations with the PT, typically ten raw bias frames 
are obtained per night. 

Once the master bias is created, master dome flats for 
each filter are made; this is done by subtracting the master 
bias from each raw dome flat frame, removing any residual 
median offset from the overscan regions, trimming off the 
MTP IPE is a suite of code written in a combination of the overscan region, applying any linearity corrections, and me- 
Tcl and C programming languages and based upon the SDSS| dian filtering the thus-processed individual dome flats for a 
DERVI SH+ASTROTOOLS (Stoughton 1995; Sergey et al. given filter. For the PT, typically five raw dome flats are ob- 
1996) software environment. The current version of MTP IPE| tained for each of the five filters (u'g'r'i'z') every afternoon 
is v8 . 3, although changes since v8 . have been mostly of scheduled observing. The same processing steps are em- 
cosmetic and have generally dealt with the smooth running ployed to create master twilight flats from the raw twilight 
of the PT Factory (§0). MTP IPE includes four main pack- flat frames. For the PT, raw twilights are obtained mostly in 
ages used in normalreductions: preMtFrames, mtFrames,|i', since we find that dome flats sufficiently flatten the PT 
excal, and kali, and they are run in that order (Fig.0. g'r'i'z' frames, whereas only twilight flats are capable of 

adequately flattening the PT u' frames. Typically 5-10 raw 
twilights (mostly in u') are obtained each evening twilight 
in which the skies are reasonably clear. 
The firstpackage, preMtFrames, is basically a "pre -burner" | Master fringe frames are needed to remove the additive 



Fig. 3 A flowchart of MT P I P E . 



and its set of SDSS filters. The SDSS magnitudes, however, 
are defined in the natural system of the SDSS 2.5m imag- 
ing camera and its SDSS filters. These two systems are very 
similar and the coefficients of the transformation equations 
are quite small. We will describe the differences in these 
two systems and provide transformation equations in §|2] 



2 The MTP IPE packages 



2.1 preMtFrames 



it creates the directory structure for the reduction of a night's 
data, including various parameter files needed as input for 
the other three packages, and it runs quality assurance tests 
on the raw data. Furthermore, it verifies information in the 
electronically generated observing log (which is based upon 
the information within the FITS image headers), identifies 
the type of each image (e.g., bias frame, dome flat, twilight 
flat, primary standard field, survey secondary patch field, or 
specially targetted manual fields), matches the frame iden- 
tifier to a list of approved standard fields, verifies that a 



effects of internal CCD illumination patterns seen in the 
i' and z' filters. To create master i' and z' fringe frames, 
mtFrames seeks out all the secondary patches observed 
during a night. These secondary patches are useful for cre- 
ating master fringe frames because they are relatively deep 
(i.e., have relatively long exposure times), and because they 
are of relatively uncrowded star fields containing few if any 
bright (highly saturated) stars, and because a given secondary! 



5 MJD is the Modified Julian Date, defined by the relation MJD : 
- 2,400,000.5, where JD is the Julian Date. 



: JD 
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Fig. 4 Histograms of pixel values for each of the raw PT bias frames taken on MJD 53501 . Note that the PT's CCD has 
2-amplifier readout electronics; hence, the bimodality of the distributions. 



patch is typically not repeat-targetted during a given night. 
If there are a sufficient number of secondary patches ob- 
served on that night (for the PT, the default minimum is 
7), master fringe frames will be created for the i' and z' fil- 
ters. The i' and z' frames for the secondary patches are fully 
processed up through flatfielding — they are bias-subtracted 
using the master bias, their residual overscan offset is sub- 
tracted, their overscan is trimmed, they are linearity cor- 
rected, and they are flatfielded with the appropriate master 
flat frame. After this processing, the secondary patch frames 
for a given filter are scaled to a default median background 
sky and median filtered. This final median-filtered image is 
the fringe frame for that filter. 

Note that, in the above description of creating master 
biases, flats, and fringes, we only make mention of data ob- 
tained over a single afternoon+night. Why not use the data 
obtained over a full week or month or more to beat down 
statistics? Actually, it is possible to do this with mtFrames,| 
if one collects all the data from a multiple-night run into a 
single directory and manipulates the ASCII observing log 
file appropriately; this process is currently not automated, 
and hence requires a certain level of human interaction. With| 
general PT operations for the SDSS, however, the guiding 
principle is that of considering each night as a self-contained 



experiment. Often, obtaining photometric calibration data 
for the SDSS is time critical, and it is ill-advised to wait 
until the end of a dark run to obtain (hopefully) all the nec- 
essary bias, flat, and fringe frames. That said, those nights 
when it is impossible to obtain a sufficient set of biases, 
flats, or (which is more often the case) secondary patches for 
fringe frames, the necessary master frames are copied from 
another night's processed data. This is permissible since the 
PT's biases and flats are stable over the course of a typical 
3-week observing run — significant discontinuous changes 
in these calibration frames are, however, often noticeable 
from one observing run to the next — and the PT i' and z' 
band fringe patterns — being an artifact of physical varia- 
tions in the thickness of the CCD itself — are stable over 
even longer periods. 

The second part of mt Frames applies these master bi- 
ases, flats, and fringe frames to the target frames — i.e., to 
the frames of fields containing primary standard stars, of 
secondary patch fields, and of any manually targetted fields 
of interest. The standard processing steps apply: for a tar- 
get frame in a given filter, the master bias is subtracted, 
any residual median offset from the overscan regions is sub- 
tracted, the frame is trimmed of the overscan regions, a lin- 
earity correction is applied, and the appropriate master flat 
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Ro* z-bona Image for Target 1406D4 Master Bias from MJD53501 




Corrected z'-band Image for Target 1406D4 




500 1000 1500 2000 



Fig. 5 The image processing steps, using PT data from MJD 53501 as an example: (a) the raw z'-band image of the 
secondary patch 1406D4, (b) the master bias frame from MJD 53501, (c) the master z'-band dome flat from MJD 53501, (d) 
master z'-band fringe frame from MJD 53466 (used for processing the PT data from MJD 53501), (e) the bias-subtracted, 
trimmed, flat-fielded, and fringe-corrected z'-band image of the secondary patch 1406D4. The (linear) stretch used in these 
images is from 90% to 1 10% of the given image's median pixel value. 
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Pre— fringe-corrected z'-band Image For Target 1406D4 Post-fringe-corrected z'-band Image for Target 1406D4 

(Central 512x512 Pixels) (Central 512x512 Pixels) 




BOO 1000 1200 800 1000 1200 



Fig. 6 An example of de-fringing a PT z'-band image: (a) the central 512 x 512 pixels of the bias-subtracted, trimmed, 
flat- fielded, but not fringe-corrected z'-band image of the field featured in Figure |5] (b) the same as (a), but after fringe- 
correction. The (linear) stretch in these two images is the same (from 95% to 105% of the median pixel value of (a)). The 
typical trough-to-peak amplitudes of the fringes in (a) is «2%. 



frame is used; if the target frame is in the i' or z' filter, a 
fringe frame is subtracted and the median sky value of pre- 
defringed frame is added back to all the pixels. 

Once a target frame is fully processed, objects are de- 
tected as peaks na above the sky background, where a is 
the rms scatter in the sky background, and aperture photom- 
etry is performed (for the general PT operations, n = 10). 
For the aperture photometry, one of two apertures may be 
used. The larger aperture is used for fields containing stan- 
dard stars and for bright stars in secondary patches and man- 
ual target-of-opportunity fields. The smaller aperture is used 
for faint stars in the secondary patches and manual fields; 
an aperture correction is applied to convert the small aper- 
ture counts into a measure of the large aperture counts for 
these faint stars. The actual sizes of these apertures are pa- 
rameters that can be adjusted. For example, when setting up 
the original standard star network (Smith et al. 2002), the 
primary standard stars were extracted using a 24-arcsec di- 
ameter aperture. This size was selected to avoid problems 
associated with defocussing the brightest stars, required for 
some of the USNO observations. 

Typically a given field is observed in multiple filters be- 
fore moving on to the next field, yielding a sequence of 
target frames for that field, one for each filter. In general 
PT operations, each target field is observed in all five filters 
(u'g'r'i' z') before moving on to the next target field. Once 
mt Frames has performed aperture photometry on all the 
frames for a given u'g'r'i' z' sequence, objects detected in 
each filter are matched (by default, they are matched to the 



detected objects in the r' filter) and the merged object list 
is written to disk as a FITS binary table file. For the PT, 
which has a known radial distortion, a distortion correction 
is applied at this stage to the aperture photometry. 

As an example of the processing steps in mt Frames, 
we continue our examination of the PT data from MJD 53501 .1 
In Figure[5] we see a raw z' image for one of the secondary 
patches observed that night (Fig.[3k), the master bias for that 
night (Fig.[5j>), the master z'-band dome flat for that night, 
(Fig.|5J;), the master z'-band fringe frame for that night (ac- 
tually taken from MJD 53466, due to an insufficient num- 
ber of secondary patches observed on MJD 53501; Fig. |5Jl), 
and the final, bias-subtracted, flat-fielded, de-fringed z' cor- 
rected frame for that secondary patch (Fig.|5^). In Figure[6] 
we zoom in on the central quarter of the pre- and post- 
fringe-corrected frame to show the effects of de-fringing. 

Finally, mt Frames outputs quality assurance plots for 
the night processed. Two examples are shown in Figure 
Figure [7^ and b show, respectively, the sky brightness (in 
ADU counts per pixel) vs. time and the mean FWHM per 
image vs. time for PT data obtained on the night of MJD 
53501. 

2.3 excal 

The third package, excal, takes the output of mt Frames 
for the u'g'r'i' z' standard star fields, identifies the individ- 
ual standard stars within these fields, and — using the in- 
strumental magnitudes for these stars as input — invokes a 
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Fig. 7 A sampling of the quality assurance plots output by mt Frames for the PT data taken on the night of MJD 53501: 
(a) Sky counts (in ADU/pixel) vs. time in each of the five passbands; (b) Mean image FWHM (in pixels) vs. time in each 
of the five passbands. (Note that the scale of the PT's CCD is 1.15"/pixel.) 
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Fig. 8 A sampling of the quality assurance plots output by excal for the PT data taken on the night of MJD 53501: (a) 
first-order extinctions (k coefficients) vs. time for each of the five passbands, and (b) magnitude residuals for the u'g'r'i'z' 
standard stars observed over the course of the night (the zero lines are staggered by 0.1 mag offsets in order to include the 
residuals from all five filters on the same plot). Note that, in (b), solid symbols denote the observations that were included 
in the solution, whereas open symbols denote the observations that were removed from the solution. 
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Extinction coefficients from photometric solution for 51182 
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Fig.9 The same excal quality assurance plots as in Figure[8] but for the USNO-l.Om data taken on the night of MJD 
51 182, in which the k terms were solved for in 3-hour time blocks. 
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Table 1 PT Color Indices (excal & kali) 



filter color index 
(x - y)o 

u ' («' - 9% 

g' (g' - r') 

r' (r'-i')o 

i' (r'-i')o 

Z (l — z') 



Table 2 PT Color Zeropoints (excal & kali) 



color zeropoint 


value 


0' - ff'kzp 


1.39 


(a' - r')o,z P 


0.53 


\T i )o,zp 


0.21 


{i Z )o,zp 


0.09 



Table 3 Instrumental Color Terms For the PT (excal) 



<MJD> 


b(u') 


Hg') 


b(r') 


b(i') 


b(z') 




Old u'g'r'i' 


z' filters 






51551 


0.001 


0.023 


0.024 


0.028 


0.002 


51580 


0.001 


0.025 


0.025 


0.029 


0.002 


51609 


0.001 


0.024 


0.024 


0.028 


0.002 


51639 


0.001 


0.023 


0.024 


0.028 


0.002 


51668 


0.001 


0.015 


0.022 


0.025 


0.002 


51698 


0.001 


-0.003 


0.017 


0.017 


0.002 


51727 


0.001 


-0.015 


0.014 


0.012 


0.002 


51786 


0.001 


-0.015 


0.014 


0.012 


0.002 


51815 


0.001 


0.008 


0.020 


0.021 


0.002 


51845 


0.001 


0.020 


0.023 


0.027 


0.002 


51874 


0.001 


0.027 


0.025 


0.029 


0.002 


51904 


0.001 


0.032 


0.026 


0.031 


0.002 


51934 


0.001 


0.034 


0.027 


0.032 


0.002 


51964 


0.001 


0.033 


0.027 


0.032 


0.002 


51993 


0.001 


0.031 


0.026 


0.031 


0.002 


52023 


0.001 


0.024 


0.024 


0.028 


0.002 


52052 


0.001 


0.014 


0.022 


0.024 


0.002 


52082 


0.001 


-0.004 


0.017 


0.017 


0.002 


52111 


0.001 


-0.016 


0.014 


0.012 


0.002 




New 


u'g'r'i 


V filters 






>52140 


0.001 


-0.041 


0.009 


0.010 


0.002 



least squares routine to calculate the photometric zeropoint 
and the atmospheric extinction in each filter passband. Since 
the overall system response in each filter (which is mea- 
sured by the photometric zeropoints) is unlikely to change 
significantly over the course of single night, a single pho- 
tometric zeropoint is calculated for a night. On the other 
hand, changes in the atmospheric transparency in each filter 
(which is measured by the atmospheric extinction term), can 
vary significantly over the course of a single night. There- 
fore, excal has the option to solve for the extinction in 
blocks of time that cover a night. In processing the origi- 
nal u'g'r'i'z' standard star network, typically 3-hour blocks 
were used; in standard PT operations, typically the block 
size is set to 15 hours or more, in order to solve for only 
a single nightly extinction in each passband. Instrumental 
color terms and second-order (color xairmass) extinctions 

www.an-journal.org 



Table 4 PT Standard Star Color Ranges (excal & 
kali) 



color 


0.70 < (u 


-a 1 


)o < 2.70 


0.15 < (<?' 


-r' 


)o < 1.20 


-0.10 < (r 


•t 

— i 


)o < 0.60 


-0.20 < (V 


f 

— z 


)o < 0.40 



may also solved for, although generally multiple nights of 
data are needed to determine these with any confidence, so 
their values are usually set to pre-determined defaults. Un- 
calibrated candidate standard stars (i.e., stars which are be- 
lieved to be non-variable but do not have previously deter- 
mined u'g'r'i'z' magnitudes) can also be used in the least 
squares solution as extinction standards; a useful output of 
the least squares routine is an estimate of their calibrated 
u'g'r'i'z' magnitudes. 

The photometric equations solved for by excal have, 
for a given filter x and color index x — y, the following 
generic form: 

X'mst *£o ~T" Ojx k x X 

+b x [(x - y) -{x- y)o,z P ] 

+c x [(x - y) Q -{x- y)o tZp ][X - X zp ], (1) 

where Xi nst is the measured instrumental magnitude in filter 
the extra-atmospheric magnitude, (x — y) is the 
extra-atmospheric color, a x is the nightly zero point, k x is 
the first order extinction coefficient, b x is the system trans- 
form coefficient, c x is the second order (color) extinction 
coefficient, and X is the airmass of the observation. X zp 
and (x — y) ,zp are zeropoint constants for the airmass X 
and the color index (x — y), respectively. 

For standard PT reductions, the filter bands are linked 
to the color indices shown in Table the color zeropoint 
constants are set to the values shown in Table[2] and X zp is 
set to 1.3 (roughly the average airmass of the standard star 
observations). Furthermore, since their values are all very 
small and since removing them simplifies the fits to the data, 
the second-order extinction coefficients (the c x 's) are all set 
to zero for standard PT reductions. 

Note that the choice of which color index is linked with 
which filter band is an option in excal. For instance, in 
our effort to extend the u'g'r'i'z' standard star network into 
the southern hemisphere (see §[5}, we link the i' band to the 
(i' — z') color index rather than to the (r' — i') color in- 
dex used in the standard PT reductions, in order to maintain 
consistency with the photometric equations used in estab- 
lishing the original u'g'r'i'z' standard star system (Smith et 
al. 2002). 

The photometric equations are solved iteratively, per- 
forming the following loop: 

1. The photometric equations are fed the current estimates 
for the extra-atmospheric magnitude, colors, and photo- 
metric coefficients (e.g., for the g' passband, these would| 
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be g' D , (g' — r') Q , a g , b g , c g , and kg). For the first iter- 
ation of this loop, initial estimates for these parameters 
are fed into the equations. 

2. The photometric equation for each of the passbands is 
solved in turn. This is done for all 5 passbands. If any 
stars (e.g., extinction standards) are having their magni- 
tudes (e.g., g' Q ) solved for, their colors (e.g., (g' — r') ) 
are kept fixed at the values they had at the first step of 
this loop. (Having the colors change as their magnitudes 
in each filter are updated as one goes from one filter's 
equation to the next has been found to yield unstable 
solutions.) 

3. Next, after the equations for all 5 passbands have been 
solved, if any extinction star magnitudes, are being solve< 
for, their colors must be updated and saved for the next 
iteration. 

4. If excal is being run interactively (a rare occurrence 
nowadays for PT processing; see §0), then the residu- 
als from the photometric solution are displayed to the 
screen, in each passband, and the user is allowed to dis- 
card stars from the solution (or add back in stars previ- 
ously discarded). 

If the user does add/remove any star to/from the solu- 
tion, excal returns immediately back to the top of the 
loop and starts again; otherwise, excal exits out of the 
loop. 

5. If excal is being run non-interactively, it sigma-clips 
the outliers in the residuals and returns to the top of the 
loop. In this case, the number of iterations of this loop 
is controlled by an adjustable parameter. For normal PT 
operations, the clipping is done at the 2.5a level, and the 
loop performs 3 iterations. 

One can also choose to fix any of the photometric coeffi- 
cients to a preset value and not solve for it. For instance, for 
the PT, we set the instrumental color terms to the preset val- 
ues shown in Table [3] (linearly interpolating between MJDs 
as needed). Note that, up through July 2001, the PT had a 
set of u'g'r'i'z filters which showed a fairly strong seasonal 
variations; the current set of PT u'g'r'i'z' filters, installed 
in August 2001, have much more stable transmissions. 

For PT reductions, where the prime focus is providing 
secondary patch stars calibrated on the SDSS 2.5m ugriz 
system, the set of standard stars used in the fits to the photo- 
metric equations are constrained to a set of relatively narrow 
color ranges (shown in Table[4}. Only standard stars that fall 
within all four of these color ranges are used, for, within 
these narrow ranges, the transformation from the USNO 
u'g'r'i'z' system to the SDSS ugriz system is best deter- 
mined and known to be very linear. 

Finally, excal, like mtFrames, also outputs quality 
assurance plots for the night processed. Figure[8]shows qual-| 
ity plots for a night of PT data (MJD 53501) in which the 
first-order extinctions (fc-terms) were solved for in one long 
(15-hour) time block that covers the whole night. Figure|5] 
shows quality plots for a night of USNO- 1.0m data (MJD 



51182) in which the k first-order extinctions where solved 
for in shorter (3-hour) time blocks. 

Figures [8^ and [9^ show the first-order extinctions for 
the two nights in question; Figures and [9j>, the resid- 
uals all of the standard stars observed over the course of 
each of these two nights. For PT reductions, we consider a 
night photometric if the rms residuals are a u > < 0.040 mag, 
o g i < 0.025 mag, a r i < 0.025 mag, oy < 0.025 mag, and 
oy < 0.040 mag. 

2.4 kali 

The last package used in normal operations, kali, per- 
forms the astrometric calibration of the secondary and ap- 
plies the photometric solutions produced by excal to the 
secondary patch instrumental magnitudes. As noted above, 
these secondary patch fields, which are observed by the PT, 
are later scanned over by the imaging camera on the SDSS 
2.5m telescope. Thus, they act as calibration transfer fields 
to set the photometric zeropoints for the imaging camera 
data. The kali package can also be applied to calibrate 
manual, or target-of-opportunity, fields observed during a 
night. 

Astrometry is performed using a triangle matching tech- 
nique based upon the one used by the Faint-Object Classfi- 
cation and Analysis System (FOCAS; Valdes et al. 1995). 6 
For standard PT processing, finding charts extracted from 
the Guide Star Catalog 1.1 (Lasker et al. 1990) are used 
for the matching. For telescopes with smaller fields-of-view 
(e.g., the USNO-l.Om and the CTIO-0.9m; see §0, how- 
ever, the deeper (and hence more densely populated) Guide 
Star Catlogue 2.2 7 is sometimes necessary to obtain a good 
astrometric solution for a given field. 

kali outputs the photometry for the secondary patches 
both in the USNO u'g'r'i'z' system and in the SDSS ugriz 
magnitudes, kali first applies the u'g'r'i'z' photometric 
solutions from excal to the instrumental magnitudes of 
the secondary patches. This step is performed iteratively, 
until the calibrated u'g'r'i'z' magnitudes have converged. 
To convert from the USNO u'g'r'i'z system to the SDSS 
2.5m telescope's ugriz system, kali applies the following 
transformation equations: 



u 


= < 




(2) 


9o 


= g' a + 0.060 x [(«/ - r') a ■ 


- 0.53] 


(3) 


r a 


= r' Q + 0.035 x [(/ - ^ , )o - 


-0.21] 


(4) 


io 


= i' Q + 0.041 x [(r'-i')o - 


- 0.21] 


(5) 




= z' Q - 0.030 x [(%' - z') a - 


-0.09]. 


(6) 




Strictly speaking, the SDSS 


ugriz maj 


>nitudes output 



by kali only apply within the narrow color ranges de- 
scribed by Table |4] although they are generally accurate for 
a much broader range, at least for normal stars (e.g., stars 

6 A standalone ANSI C version of this code can be obtained from 

http :// spiff . rit . edu /match/ 

7 |http : //www-gsss . stsci . edu/gsc/gsc2/GSC2home.htrJ 
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Percentage of stars matched in MT secondary fields for run 53501 
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Fig. 10 One of the quality assurance plots output by kali for the PT data taken on the night of MJD 53501: the 
percentage of stars detected in the r'-band matched to stars in the finding chart file for each of the secondary patches and 
manual target fields observed during the night (note that the finding chart file for each secondary patch and manual target 
field is extracted from the the Guide Star Catalog-I; Lasker et al. 1990). 



showing no strong emission features in the their spectra) 
blueward of the M0 spectral type (see, for example, Fig. 7 
and8of Rideretal. 2004). 

Further, like mt Frames and excal, kali outputs some| 
quality assurance plots for the night processed, a sampling 
of which is shown in Figures ^| and ^2 Figure [Tol shows 
how well the PT patches were matched to the the stars in 
the finding charts, and Figures ITTk. b, and c show the ugr, 
gri, and irz color-color diagrams, respectively, for the stars 
in the secondary patches observed on this night (the ugr and 
gri diagrams also show a stellar locus (Lenz et al. 1998) for 
comparison). 

Finally, kali automatically tags the quality of the pho- 
tometry within individual secondary and manual patches by 
monitoring the position of certain features in each patch's 
color-color diagrams. For determining the quality of the pho-| 
tometry in the griz bands, kali uses the location in g 1 — r 1 , 
r' — i', i' — z' space where the red and blue branches of the 
stellar locus cross — i.e., roughly where stars of spectral 
type M0 sit, which is at g' —r' w 1.35, r' — i' « 0.50, and 



i' — z' w 0.25; this is the so-called Zhed point test, a more 
refined version of which is described in detail by Ivezic et 
al. (2004). The quality of the u band photometry is automat- 
ically tagged based upon the estimated u' — g' color of stars 
in the blue branch at g' — r' = 0.50, which is v! — g' w 1.45. 

2.5 Auxiliary Packages 

In addition to the above four packages, which are used for 
normal night-to-night survey operations, there are two other 
packages — solve_network and superExcal — which| 
were developed and used for setting up the u'g'r'i'z' pri- 
mary standard star network of Smith et al. (2002). These 
two packages are similar to excal, in that they invoke 
least squares routines to solve for the photometric param- 
eters of a set of data and for the best-fit magnitudes of the 
primary standards . They differ from e x c a 1 in that, whereas 
the main purpose of excal is the standard, single-night 
solution of the photometric equations during routine sur- 
vey operations, the main purpose of solve_network and 
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Fig. 11 Color-color quality assurance plots output by kali for the PT data taken on the night of MJD 53501: (a) the 
combined g — r vs. u — g color-color diagram, (b) the combined r — « vs. g — r color-color diagram, and (c) the combined 
i - z vs. r - i color-color diagram for the stars in all the secondary patch and manual fields from the night. The lines in (a) 
and (b) show the u'g'r'i'z' stellar loci (based on the stellar loci of Lenz et al. 1998). 
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superExcal was to calibrate the SDSS standard star net- 
work. Both solve_network and superExcal use a sin-( 
gle star to set the zeropoint for the photometric solution. 
For setting up the SDSS standard star network, this star was 
the F subdwarf BD+17°4708, whose magnitudes were de- 
fined to be u' = 10.56, g' = 9.64, r' = 9.35. %' = 9.25, 
z' = 9.23, based upon the spectrophotometric calculations 
of Fukugita et al. (1996). Thus, BD+17°4708 sets the zero- 
point for the entire standard star network. 

Although they share the same goal and the same ba- 
sic methodology, solve_network and superExcal do 
differ in some important respects, the most important be- 
ing that the code for solve_network was constructed 
independently of that for excal, whereas superExcal 
was basically an outgrowth of excal. As a result, the sep- 
arate results from these two different packages provide use- 
ful cross-checks on the final calibrations of the standard 
star system. Results and output of solve_network and 
superExcal can be found in Smith et al. (2002). 

Finally, although not strictly an auxiliary package but 
more of an auxiliary mode of operations is the MTPIPE 
follow mode, which permits the near-realtime running of 
mt Frames and excal during data acquisition at the PT. 
MTPIPE in follow mode was used by the SDSS observers 
during the 2000-2001 time frame in conjunction with the 
hoggpt photometricity monitor (Hogg et al. 2001) in or- 
der to track sky conditions over the course of a night. This 
functionality is nowadays completely handled by the dedi- 
cated hoggpt photometricity monitor. 

3 Tests of mtpipe 

As a critical component in the photometric calibration of 
the SDSS, MTPIPE has undergone a variety of tests to en- 
sure the validity of its outputs. Most of these tests have been 
aimed at the potentially trickier parts of MTP IPE — the im- 
age processing and aperture photometry in mt Frames and 
the photometric equations solver in excal. 

A major test of mtFrames was accomplished by com- 
paring its outputs with those of the hoggpt APO photo- 
metricity monitor (Hogg et al. 2001) for various nights of 
PT data. The hoggpt software was written completely in- 
dependently of mtpipe in the data analysis language IDL. 8 
It thus could act as an independent check of the data pro- 
cessing methods utilized in mtFrames. Overall, it was dis- 
covered that these pipelines created virtually identical mas- 
ter bias frames, master flat fields, and master fringes. Fur- 
thermore, for a given aperture size, the two pipelines yielded 
essentially indistinguishable aperture photometry. Plots com-| 
paring the aperture photometry from the two pipelines are 
shown in Figures [Q] El and^l In these plots, an aperture 
size of 7.43"radius is used on PT data from MJD 51639. 

The basis of the excal photometric equations solver 
is an iterative non-linear matrix inversion routine which si- 

8 |http : / /www . rsinc ■ com/| 



multaneously solves for all parameters in each filter's pho- 
tometric equation. To test this, a very stripped-down ver- 
sion of excal ("subExcal") was written which utilized 
a vanilla linear least squares line-fitting routine based upon 
the L INF IT routine of Bevington (1969) converted into the 
SDSS Tcl/ASTROTOOLS environment. subExcal does not! 
solve for all parameters simultaneously, but first solves for 
the photometric zeropoint and the first order extinction (the 
a and k coefficients) by fitting a straight line to an equation 
of the form (taking the g' filter as our example): 

Si'nst " 9o = a 9 + k 9 X ' ( ? ) 

(see Fig. I15i . and then uses the values of the a and k terms 
thus obtained to solve for the b coefficient by fitting another 
straight line to an equation of the form: 

9inst - 9o - a g - k g X = b g (g ' - r') (8) 
(see Fig. I16i. In Table|5] we compare the output of excal 
against the output of subExcal for PT data from the night 
of MJD 51639. Note that in many cases, the results are in- 
distinguishable, and in all cases the values differ by no more 
than 1(7. 

The kali package of MTP IPE is more or less straight- 
forward. It performs astrometry on secondary patches, and 
the astrometry needs to be accurate enough for PT secondary! 
patch stars to be cross-matched with stars in the SDSS 2.5m 
imaging camera data. Residuals in the astrometric solutions 
(based against the Guide Star Catalog-I; Lasker et al. 1990) 
are typically less than 1" RMS in both RAxcos(DEC) and 
DEC, which is sufficient for the PT-2.5m cross-matches. 
Further, kali applies the photometric solutions from excal| 
to the instrumental magnitudes of the secondary patches. 
This has been verified by manually applying the excal so- 
lutions to the instrumental magnitudes of the PT secondary 
patch data and comparing the results; the resulting differ- 
ences are typically of the order of the roundoff error of the 
calculator/computer used for the manual computation. 

Of course, in the final analysis, the true-end-to-end test 
of MTP IPE is the resulting quality of the SDSS photometric 
calibrations. As stated in the introduction of this paper, the 
rms errors measured over the several thousand square de- 
grees that currently make up the SDSS imaging data set are 
«0.02 mag (2%) rms over several thousand square degrees 
of sky in g, r, and i, and «0.03 mag (3%) in u and z (Ivezic 
et al. 2004; Adelman-McCarthy et al. 2006). 

In the future (§0, MTP IPE will likely need to calibrate 
secondary patches at low Galactic latitudes for SDSS-II. 
How close to the Galactic Plane (i.e, at what stellar densi- 
ties) can MTPIPE — which does aperture, not PSF, pho- 
tometry — accurately calibrate these patches? One of us 
(CTR) is currently testing MTP IPE against DAOPHOT (Stet- 
son 1987) for some low-Galactic-latitude PT fields. 

4 The PT Factory 

For the photometric calibration of the SDSS 2.5m imaging 
camera data, a method is needed for the timely, end-to-end 
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Counts(hoggpt)/Counts(mtpipe) vs. Time for apo20— mjd51 639 
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Counts(hoggpt)/Courits(mtpipe) vs. USN040 Magnitude for apo20-mjd51 639 
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Counts(hoggpt)/Counts(mtpipe) vs. USN040 Color for apo20— mjd51 639 
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"subExcal": mtpipe instrmag — stdmag vs. airmass for mjd51639 
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subExcal": mtpipe (instrmag — stdmag - a — kX — c*(color— color0)*(X-X0)) vs. color for mjd51639 
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Table 5 Comparison of excal and "subExcal" for 
MJD 51639 



Filler 


Method 


a filt.nr 


^filtr^r 


*filt.rr 




excal 
subExcal 


-19.286 ± 
-19.318 ± 


042 
.037 


-0.011 ± 0.014 
-0.011 ± 0.013 


0.516 ± 0.025 
0.530 ± 0.025 


9' 


excal 
subExcal 


-20.668 ± 
-20.673 ± 


019 
.021 


+0.026 ± 0.007 
+0.026 ± 0.007 


0.184 ± 0.012 
0.195 ± 0.014 




excal 
subExcal 


-20.578 ± 
-20.573 ± 


.021 
022 


+ 0.033 ± 0.015 
+0.033 ± 0.015 


0.119 ± 0.013 
0.118 ± 0.014 


.1 


excal 
subExcal 


-20.194 ± 
-20.188 ± 


.020 
.020 


+0.026 ± 0.023 
+0.025 ± 0.022 


0.071 ± 0.013 
0.068 ± 0.013 




subExcal 


-19.162 ± 
-19.164 ± 


026 
.025 


-0.025 ± 0.029 
-0.025 ± 0.029 


0.067 ± 0.017 
0.067 ± 0.017 



processing of the PT data that is fully integrated and main- 
tained within the official SDSS data processing environ- 
ment at Fermilab. This need is fulfilled by the PT Factory, 
a term which applies to a suite of scripts that oversees the 
day-to-day processing and quality assurance of PT data at 
Fermilab. These scripts reside in the dp software data pro- 
cessing product, which is a collection of UNIX shell scripts 
and Tcl/ASTROTOOLS code that oversees the overall pro- 
cessing and quality assurance for 2.5m imaging, 2.5m spec- 
troscopy, and PT data (Stoughton et al. 2002b). 

Each morning, the PT Factory performs an r sync 9 copy| 
of the previous night's PT data from APO to a local disk 
at Fermilab (typically ^2 GB for a photometric night). The 
morning r sync also copies the night's images and log from 
the APO 10 micron All-Sky Camera 10 (Hogg et al. 2001); 
the log, which is a record of the sky rms for each all-sky 
camera image over the course of a night, is used by the PT 
Factory to flag and exclude from mtpipe processing any 
PT target images afflicted by clouds. 

Once the PT Factory determines that the night's PT data 
have been completely transferred (typically by the early- 
to mid-afternoon of the same day), MTPIPE is automati- 
cally run on the data. Once MTP IPE has run to completion, 
there is a manual inspection step, in which a human checks 
over the quality assurance plots generated by mtFrames, 
excal, and kali. This manual inspection step exists to 
ensure that nothing unexpected occurred with either the raw 
data or with the data processing that was not already caught 
by MTPIPE's automated quality assurance tests. If, as is 
generally the case, the night's processing passes the man- 
ual inspection step, the photometric solutions and survey- 
quality secondary patches are copied to a location on disk 
where NFCALIB expects to find them, and the contents of 
the night's reductions are loaded into a database and saved 
to a tape robot. 

An additional input into the PT factory is the "ptlog" 
observer's report, which is e-mailed at the end of each night 
to a ptlog mailing list. This report contains summary com- 
ments on the night, including the sky conditions, which (if 
any) manual or special targets were observed, any issues 
with the PT hardware or observing software, and any other 



9 |http : / / r sync . samba ■ org/| 

10 |http : //hoggpt ■ apo ■ nmsu ■ edu/irsc/ tonight /| 



comments that might be relevant to either the PT data pro- 
cessing or to the PT operations. The ptlog is read each morn- 
ing, and its contents are used in deciding whether any man- 
ual intervention of the PT Factory is needed for that night's 
data. As such, it provides an essential link between the ob- 
servers at APO and the data processors at FNAL. 

A final step of PT processing, not strictly part of the 
PT Factory, is updating the secondary patch database main- 
tained at APO. This database keeps track of which secondary! 
patches have been observed, which have not been observed, 
and which of the observed patches have been determined 
by the PT Factory to be of survey quality. It also keeps 
track of the observing priority of the secondary patches. 
This database therefore serves as essential feedback to the 
observers regarding the secondary patch observing priori- 
ties for a given dark run or even a given night. 

5 Other SDSS-related Projects using 

MTPIPE 

Since the inception of the SDSS observing effort on the 
PT, we were aware that MTPIPE would enable additional 
projects to be undertaken. Initially, these were centered on 
use of the PT during bright time, but later it became appar- 
ent that other telescopes could be used with modification of 
the preMtFrames code and creation of additional input 
parameter files. 

Early in the SDSS we recognized the usefulness of open 
star clusters as potential calibration fields for the main sur- 
vey, as well as being good science targets in their own right. 
This led to the development of a star cluster project in the 
u'g'r'i'z' filter system. Also, shortly after the establishment 
of the original u'g'r'i'z' standard star network (Smith et 
al. 2002), we began to receive requests for southern hemi- 
sphere standard stars. These requests led us to develop a 
southern u'g'r'i'z' standard stars project which was under- 
taken as an NOAO Survey Program. Further, as we reached 
the point where we could see the end of the SDSS proper, 
development of survey extension strategies were pursued by 
the collaboration. The PT and the MTPIPE software are in- 
tegral to the support of these extension efforts. 

In the rest of this section, we will briefly discuss each of 
the project areas that are enabled by this unique calibration 
software, and point to references for additional information. 

5.1 Southern u'g'r'i'z' Standards Project 

Due to the nature and location of the SDSS, the bulk of 
the initial set of standard stars to calibrate the survey were 
placed in the northern hemisphere and on the equator (Fig.Ql.| 
As part of the NOAO Surveys Program, 1 1 a group of us was 
granted time on the CTIO-0.9m over the 2000-2004 time 
period in order to establish a network of u'g'r'i'z' standard 
stars in the southern hemisphere (see Fig.Q. Processing and 



11 |http : / / www . noao ■ edu/ gateway / surveys /programs ■ htm J 
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analysis of these data are nearly complete, and we are tar- 
getting the second-half of 2006 for release of the bulk of the 
data. (An initial data release for standard stars in the Chan- 
dra Deep Field South appeared in Smith et al. 2003.) As the 
processing and analysis is completed, results will be made 
available on our publicly accessible website. 12 

5.2 u'g'r'i'z' Open Clusters Project 

One of the most powerful tools available in astrophysics for 
exploring and testing theories of star formation, stellar and 
galactic evolution, and the chemical enrichment history of 
the Galaxy is the study of open star clusters. Therefore, 
a group of us embarked on a survey of (mostly) southern 
hemisphere open star clusters using the u'g'r'i'z' filter sys- 
tem on small telescopes. This project had as its original 
goal to verify the SDSS calibration scheme. This goal has 
since evolved. Current plans for observing star clusters in 
the SDSS filter system now include: (1) their potential use 
in calibrating the SDSS extension programs, and (2) creat- 
ing a large uniform imaging survey of clusters similar to the 
WIYN Open Cluster Survey (WOCS) (see Mathieu 2000). 

In the first of our cluster papers (Rider et al. 2004), we 
present results for NGC 2548. We chose this particular clus- 
ter because we had observations of it using three telescopes: 
the CTIO Curtis-Schmidt, the PT, and, most importantly, the 
USNO-1 m, the telescope with defines the u'g'r'i'z' pho- 
tometric system (Smith et al. 2002). A second cluster pa- 
per (Moore et al. 2006) has been submitted and addresses 
NGC 6134 and Hogg-19. Further papers are in preparation. 
In all, we have data in hand for «100 open clusters. 

These data may be used to verify the recent age and 
metallicity models of Girardi et al. (2004) and the prior 
work of Lenz et al. (1998), and to verify and expand upon 
the u'g'r'i'z' to UBVR C I C transformations presented in the 
standard star paper (Smith et al. 2002). Rodgers et al. (2006) 
is taking another look at these filter system transformations 
for main sequence stars in an effort to sort out luminosity 
class effects. 

5.3 SDSS-II 



order to study the structure and evolution of the Milky 
Way (Newberg et al. 2004). 
3. The SDSS-II Supernova Survey, which will image 250 
sq deg along the Celestial Equator every two nights dur- 
ing the Autumns of 2005-2007, with the aim of obtain- 
ing well-sampled, well-calibrated, multi-band lightcurves| 
for 200 Type la supernovae in the redshift "desert" z = 
0.1 - 0.35 (Frieman et al. 2004; Sako et al. 2005). 

The PT and MTPIPE are playing an important role in all 
three of these components of the SDSS-II. First of all, since 
it is basically a continuation of the original SDSS within the 
original survey borders, the Legacy component of SDSS- 
II continues to benefit from the use of MTPIPE reductions 
of PT secondary patches as part of the general photometric 
calibration strategy of the 2.5m imaging camera data. 

Likewise, the 2.5m imaging scans for SEGUE also uses 
MTP IPE reductions of PT secondary patches located within 
the SEGUE survey area (see Fig.|2jl. But the PT and other 
small telescopes (like the USNO-l.Om) have an even big- 
ger role to play in SEGUE. One of the major thrusts of 
SEGUE is the SEGUE Open Cluster Survey. This is a wide- 
field imaging program with limited spectroscopic follow-up 
within SDSS-II. While the imaging program will not go as 
faint as the WOCS (see Mathieu 2000), we will observe al- 
most an order of magnitude more clusters than the WOCS. 
Directed pointings for specific calibration targets and for 
the brighter stars of the clusters will be undertaken with 
the smaller telescopes, and those data will be reduced with 
MTPIPE. 

Finally, the PT and MTPIPE have a role to play in the 
Supernova component of the SDSS-II. In this case, the role 
has less to do with photometric calibrations — since the pri- 
mary Supernova strategy is to use differential photometry of 
repeat 2.5m imaging scans against previously calibrated ar- 
eas of the sky scanned by the original SDSS survey — but 
more of a role in followup observations. Here, the PT will 
be able to follow the light curves of the brighter candidate 
Type la supernovae detected by the main SDSS-II Super- 
nova Survey during those times when 2.5m telescope time 
is allotted either to the Legacy survey or to SEGUE, and 
MTPIPE will be able to process these PT observations. 



The SDSS officially ended on 2005 June 30. The SDSS-II 
is a 3-year extension of the SDSS beyond this date, and is 
composed of three parts: 

1. The Legacy Survey, which will fill in a gap in the SDSS 
imaging and spectroscopic coverage of the Northern Galac 
tic Cap due to worse-than-expected weather conditions 
during the five years of the original SDSS. 

2. The Sloan Extension for Galactic Understanding and 
Exploration, or SEGUE, which will image about 3500 
sq deg of sky — mostly at Galactic latitudes of \b n \ < 
30° — and measure spectra to determine abundances 
and radial velocities of a quarter of a million stars, in 

12 |http : / /www-star . f nal ■ gov71 



6 The Future 

MTPIPE is a mature and stable data processing pipeline. 
The last major revision — from v7 . 4 to v8 . occurred 
(ver four years ago (on 2002 April 22). As noted in § [2] 
the current version is v8 . 3, and the changes since v8 . 
have been mostly cosmetic and have generally dealt with 
the smooth running of the PT Factory (§0}. This being the 
case, what is the future of MTP IPE? 

There are, of course, many relatively small changes to 
MTPIPE that could yield incremental improvements in the 
final output and to the photometric calibrations of the SDSS. 
These include updates to the de-fringing code in mtFrames| 
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and continued improvements to the automated quality as- 
sessment of secondary patches for the PT Factory. Also, 
improving the compatibility of the output of MTPIPE with 
other astronomical data packages would be a plus. Along 
these lines, adding World Coordinate System (WCS; Cal- 
abretta & Greisen 2002; Greisen & Calibretta 2002) key- 
words to the FITS headers of MTPIPE corrected images is 
high on the list of priorities. (Currently, there is a standalone 
Tcl/ASTROTOOLS script to do this, but it has yet to be fully 
incorporated into MTP IPE.) 

A larger task that may be ahead for MTPIPE is how 
to deal with crowded field photometry, mt Frames is an 
aperture photometry code, but many of the science projects 
that currently use or plan to use MTPIPE involve crowded 
fields — in particular, star clusters (§|5). Furthermore, the 
photometric calibration SDSS-II/SEGUE scans may require 
the placement of SEGUE secondary patches at very low 
Galactic latitude. With the task of crowded field photom- 
etry in mind, then, we are investigating the possibility of 
adding a toggle switch to mt Frames which would allow 
one to choose to run a PSF-fitting photometry package like 
DAOPHOT (Stetson 1987) or DoPHOT (Schechter et al. 
1993) on corrected frames output by mt Frames; a new 
script, which would be run after mt Frames, would then 
package the resulting DAOPHOT/DoPHOT output into an 
excal-/kali-ingestible format. 
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